clear all;
clc;
% Run the Monte Carlo simulation on a given number of photons
% with various medium properties.


% Specify the number of photons to simulate and the various 
% properties to use during the simulation.
No_photons=1000;

mu_a = 0.1;
mu_s = 100;
n1=1;                   %Refractive Index of Air
n2=1.37;
g = 0.9;
n = 1.37;

for(i=1:1:No_photons)
    Pencil_Beam=MonteCarlo(mu_a, mu_s, n1, n2, n, g);
    plotphoton(Pencil_Beam);
    hold on;
end
